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Analyses of the Earth Radiation Budget Experiment (ERBE) data show that the 
effects of clouds on the solar and thermal infrared radiation in the tropical deep 
convective regions have a similar magnitude but opposite signs. 

This small difference in the effects of clouds on radiation led Hartmann et al. 
(2001) to conclude that the contrast in the net radiation at the top of the 
atmosphere between the convective and non-convective regions must also be 
small. However, we have found that the ERBE data do not generally show a 
small contrast in the radiation between the convective and non-convective 
regions, and the model used by Hartmann et al., therefore, seems unlikely to 
represent the real physical processes involving convection, radiation, and climate 
in an appropriate way. 
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Analyses of the Earth Radiation Budget Experiment (Barkstrom, 1984) data show that the 
shortwave (SW) and longwave (LW) cloud radiative forcing (CRF) at the top of the atmosphere 
(TOA) in the tropical deep convective regions have a similar magnitude but opposite signs. The 
CRF is defined as the difference between the full-sky and clear-sky radiative fluxes. Hartmann et 
al. (2001) (referred as HMF hereafter) computed the CRF for individual cloud types using a 
radiative transfer model with the cloud information taken from the International Satellite Cloud 
Climate Project (ISCCP) (Rossow and Schiffer, 1991) data archive. They found that individually 
each cloud type had a large impact on the radiation at the TOA, but collectively the effect was 
small. From the analyses of the ERBE data, they claimed that the net CRF at the TOA is small 
not only over deep convective regions but also over the non-convective regions. The perceived 
small CRF led them to conclude that the contrast in the net TOA radiation between the convective 
and non-convective regions must also be small. HMF used a simple model to investigate the 
relation between convection and radiation in the tropics. They hypothesized that the small 
contrast in the TOA radiation between the convective and non-convective regions was due to the 
high sensitivity of circulation to the sea surface temperature (SST) gradient and the high 
sensitivity of cloud albedo to vertical velocity. 

However, the ERBE data do not generally show a small contrast in the TOA radiation 
between the convective and non-convective regions, and the model used by HMF, therefore, 
seems unlikely to represent the real physical processes involving convection, radiation, and 
climate in an appropriate way. 

Figure la shows the net radiation and CRF in the Pacific averaged over five years (1985- 
1989) for the equatorial band 10°S-10°N derived from the ERBE data. The net heating of the 
earth-atmosphere system attains a maximum value of -80 W m ^ west of 160°E and decreases 
continuously eastward to a value of -50 W at 260“E. The longitudinal distribution of the CRF 
(Figure lb) is very similar to that of the net radiation. The CRF ranges from -5 W west of the 
dateline to -32 W at 260°E. The contrast between the deep convective region west of 160°E 
and the suppressed region to the east is 1; 1.6 for the net radiation and 1:6 for the CRF. These 
contrasts can hardly be considered as small. Averaged over the five ERBE years (1985-1989), the 
full-sky net TOA radiation in the tropics (30°S-30°N) is 48.6 W the clear-sky net radiation is 
61.1 Wm \ and the CRF is -12.5 Wm^. Thus, the CRF in the tropics as a whole is significant. 
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HMF suggested that the small CRF in the deep convective region as found by Kiehl (1994) 
is not a coincidece but is a combined effect of an ensemble of various types of clouds occurring at 
various heights and v^ith a large range of optical thickness. They further suggested that the small 
CRP is a result of the feedback between convection, cloud albedo, and the SST gradient. In fact, 
the small CRF in the deep convective region is easy to explain. Clouds have the effects of 
reflecting SW radiation and reducing OLR. Except for thin clouds at high altitudes, the SW CRF 
in the tropics is greater than the LW CRF. This can be seen in Fig. lb that the net CRF at 260 E is 
large, approximately -30 W rri^. This region is in the descending branch of the Walker circulation 
where the deep convection is suppressed and the boundary-layer clouds are extensive. For those 
boundary-layer clouds, the contrast between the cloud-top temperature and the SST is small, and 
so is the LW CRF. On the other hand, those clouds are highly reflective to SW radiation, and the 
SW CRF is large. The result is a large negative CRF. Westward from 260"E, the boundary-layer 
depth and, hence, the cloud height increase, but the cloud cover decreases. The increase in cloud 
height enhances the difference between the cloud temperature and the SST, which leads to an 
enhanced LW CRF, Thus, the net CRF decreases westward toward the dateline. West of the 


dateline is the deep convective region, where high-level clouds are extensive. The CRF is small, 
within -10 W m^. 

In the simple model of HMF, the tropics is divided into a convective region and a non- 
convective region. Each of these two regions is characterized by the SST, T, the net radiation at 
TOA, /?, and the areal extent, A. The rate of change of the SST is given in their model by 
M 


at At 


dT2 

dt 


= R2+-^cje,i-d,2)+Exp2 


( 1 ) 

( 2 ) 


where c is the heat capacity of the ocean mixed layer, M is the mass flux between the two regions, 
Cp is the heat capacity of air, 8^ is the saturation potential temperature of the surface air, Exp is the 


export of heat out of the domain of interest, and the subscripts 1 and 2 denote the convective and 
non-convective regions, respectively. 

Assuming that exports of heat out of the convective and non-convective regions are the 
same, equilibrium solutions for T;, T 2 , Ri, and R 2 are obtained by time integration of Eqs. (1) and 
(2), so that 




1 1 


( A] A2 J 


(3) 
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Thus, the imbalance betv/een Rj and R 2 equals the transport of heat from the convective region to 
the non-convective region. 

The mass flux between the two regions, which is related to the vertical pressure velocity co. 


is given by 


6)9 , CO ] 

A/ = A2 — = — 


(4) 


where g is the gravitational acceleration. The convective area fraction Aj is set to be 0.3. The 
vertical pressure velocity in the convective region is assumed to be linearly proportional to the 
temperature contrast between the convective and non-convective regions, 

(0^ = -y(Ti-T,)^-YAT ( 5 ) 


where /is a disposable constant. 


The net radiation at TOA is given by 

where S is the insolation at TOA, oc is the albedo, F is the OLR, and i=l or 2. The OLR in each 

region is assumed to vary only slightly with the SST and is practically fixed at 190 W m ^ and 280 
W m'^ in the convective and non-convective regions, respectively. The albedo in the non- 
convective region is fixed at 0.1, and the albedo in the convective region is assumed to be a 

function of vertical pressure velocity and SST, 

cCi=a2-^C0ig*{Ti) O) 

where A is a disposable constant, and q’ is the saturation water vapor mixing ratio. 

Equation (5) has the convection and, hence, atmospheric overturning increasing with 
increasing temperature contrast between the two regions. Equation (7) has the cloud albedo in the 
convective region increasing with increasing convection for a positive A. As a, increases with 

increasing T,, the net radiation R, will decrease, which will then lead to a decrease in both T, and 
the temperature contrast between the two regions. The decrease in the temperature contrast will 
cause weaker convection, which, according to Equation (7), will reduce the albedo (Xj for a 


positive A. Thus, it constitutes a negative feedback loop between SST, convection, cloud albedo, 
and radiation. The strengtli of the feedback is dependent upon the values of /and A. 


From Equation (7), we have 
1 

q\T0 Anil 


( 8 ) 
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To estimate the albedo sensitivity parameter A, HMF derived values of and AcOj by taking the 

albedo difference and the vertical velocity difference between the convective and non-convective 
regions, 

Aa, = aj-a2 

( 10 ) 


so that 


AcOi = coi -CO2 


A = - 


q(T0 


ai-a2 

co^-a>2 


( 11 ) 


By assuming an albedo contrast of 0.25 between the convective and non-convective regions, they 
found 


= -1.25x10^ (Pa 


( 12 ) 

(O1-CO2 

and A=5xl0^ (Pa s‘) ‘ for q*=0.Q2. Thus, A has a positive value, and the albedo in the convective 
region increases as convection strengthens. 

Since, as applied by HMF, and Ao), are the changes in albedo and vertical velocity of 

the convective region between either two different climate states or two different times. Equations 
(9) and (10) are of course incorrect. Instead, the correct representations should be 

Aai = aj"-a|* (13) 


Aft), = (Oi -o>| 


(14) 


where the superscripts a and b denote two different climate states. Therefore, the sensitivity of 
albedo to the vertical velocity estimated by HMF based on the differences between the convective 
and non-convective regions is unrelated to the sensitivity given in Equation (8). Equation (8) 
considers only the albedo and vertical velocity in the convective region. There is no evidence that 
albedo in the convective region increases with increasing vertical velocity as HMF claimed. 

From Equations (5) and (7), HMF estimated the sensitivity of albedo to SST from 

Aa 


AT 




(15) 


Based on the results of Ramanathan and Collins (1991) on the relation between the interannual 
variations of the SST and the absorption of solar radiation, HMF estimated that 
Aa 


AT 


= 0.04 K ' 


(16) 
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and A/-2 for ^*=0.02. They found that for Ay=2 the model predicted a strong feedback that 


brought the net radiation in the convective region close to that in the non-convective region. 

There are problems with this analysis of the albedo sensitivity to SST. Firstly, Equation 
(15) is not correct. From Equations (5) and (7), the albedo sensitivity should be 
Aai 


A(ri-72) 


^Xyq 


(17) 


so that the albedo in the convective region changes with the temperature contrast between the 
convective and non-convective regions. 

Secondly, the sensitivity of albedo to SST cannot be derived by comparing the data 
associated with interannual variations. To demonstrate this point, we show in Figure 2 a scatter- 
plot that relates albedo to SST in the central equatorial Pacific (10"S-10‘'N, I70“E-150°W) during 
April 1985 and April 1987. These two months are the same as that used in Ramanathan and 
Collins (1991). Each point represents monthly mean albedo from ERBE data archive and SST 
from the National Centers for Environmental Prediction (NCEP) analysis (Reynolds and Smith, 
1994). The spatial resolution for each point is 2.5"x2.5" latitude-longitude. The year 1987 is an El 
Nino, and the SST in the central and eastern equatorial Pacific was abnormally high. Figure 3 
shows the longitudinal distribution of equatorial SST of these two months. Compared to April 
1985, the SST in the central and eastern equatorial Pacific increased significantly by 1.0-1.4“C in 
April 1987, but the warm pool region (140“E-160“E) cooled slightly by ~0.2"C. 
Correspondingly, the center of deep convection shifted eastward by ~50‘* longitude from the 
maritime continents to the central equatorial Pacific (Chou, 1994). Associated with the eastward 
shift of the deep convecUon, Figure 2 shows that both albedo and SST in the central equatorial 
Pacific are larger in April 1987 than in April 1985. These changes are unrelated to the sensitivity 
given in Equation (17) and are irrelevant to the model of HMF. In fact, the central equatorial 
Pacific was in the descending domain of the Walker circulation during April 1985 but was in the 
ascending domain during April 1987. The albedo difference between these two months is 
equivalent to the incorrect equation (9) but not equivalent to the correct equation (13). 

Strong updrafts in tlie convection cores constitute only a small fraction of the tropics. Most 
of the tropics are in the regions where the air is descending. Even the detrained thick anvil clouds 
are in the regions of descending motion. The simple model of HMF assigns a fractional area of 
0.3 to the convective region. With this large extent of the convective region, it must include both 
the deep convection cores and the extended anvil and cirrus clouds, both optically thick and thin. 
It is not clear what the albedo of the convective region in Equation (7) refers to. If it refers to 
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the albedo of the deep convection cores, then the albedo can change very little because the cloud 
particle content and optical thickness are already high. If the convective region refers to the 
region comprising both the convection cores and the anvil clouds, then the albedo of this 
convective region might decrease with increasing SST and the gradient of SST due to the area 
response as suggested by Lindzen (1997). In Lindzen et al. (2001), the fractional cover of high- 
level clouds due to the detrainment of deep convective clouds decreases with increasing SST. 
Recently, Wielicki et al, (2002) reported that the satellite-inferred OLR in the tropics increases 
during the past decade. They suggested that this increase of OLR could not be accounted for by 
the change in dear-sky radiation and was due to a decrease in clouds. Results from those studies 
imply that the mean albedo in the convective region of the HMF model might not necessarily 
increase with increasing Ty and cOj, Furthermore, those studies also imply that the change in high- 

level cloud cover could have a large impact on the LW radiation and, hence, climate. This LW 
impact was essentially ignored by HMF. 

Finally, Equations (1) and (2) are incorrect. The terms involving the rate of SST change are 
put there for the purpose of obtaining equilibrium solutions by time integration. For equilibrium 
solutions, the /?-terms balance the M-terms as shown in Equation (3). If the mass flux M is 
defined only for the atmosphere, then Equations (1) and (2) are correct only if the net surface heat 
flux is negligible in both the convective and non-convective regions. In this case the net radiative 
heating of the earth-atmosphere system R is the same as the net heating of the atmosphere. 
Although the net surface heating has been found to be small in the Pacific warm pool during 
northern hemispheric winter (Chou, et al., 2000), it is not generally so in other tropical regions 
and seasons. Otherwise the heat fluxes at the sea surface in the tropics would not be an issue in all 
long-term climate studies. Figure 4 shows the net downward radiation at the top of the 
atmosphere, the net surface heating, and the net atmospheric heating in the region (30'’S-30‘*N, 
lOO'^E-lVO'^W) during January-August 1998. The TOA net radiation is taken from the Tropical 
Rain Measuring Mission (TRMM) Clouds and the Earth's Radiant Energy System (CERES) 
archive (Wielicki, et al., 1996). The radiation component of the net surface heating is inferred 
from Japan's Geostationaiy Meteorological Satellite radiance measurements (Chou et al., 2001). 
Whereas the surface turbulent heat fluxes are derived using the NCEP SST analysis and the 
Special Sensor Microwave/Imager (SSM/I) precipitable water and wind retrievals (Chou, et al., 
1997). The atmospheric heating is the difference between the TOA radiation and the surface 
heating. It can be seen that the net surface heating has a significant spatial variation, ranging 
approximately from -60 W to -h60 W m^. As a result, the net heating of the atmosphere has 
little resemblance to the TOA radiation. 
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On the other hand, if the mass flux M is defined as the total mass transport in both the 
atmosphere and the ocean, then to can no longer be the vertical pressure velocity of air. In this 

case. Equation (7) and aU discussions of to in HMF become irrelevant. 
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Fi gure Captions 


Figure 1, Longitudinal distributions of the net radiation (upper panel) and cloud radiative forcing 
(lower panel) in the equatorial Pacific (lO^S-lO^N) averaged over five years (1985-1989). 
The radiative fluxes are taken from the ERBE data archive. 

Figure 2. A scatter-plot relating albedo to the sea surface temperature in the central equatorial 
Pacific (10“S-10‘N, 170‘’E-150‘’W)) during April 1985 and April 1987. Each point 
represents monthly mean albedo from ERBE archive and SST from the NCEP analysis. 
The spatial resolution for each point is 2.5“x2.5‘’ latitude-longitude. 

Figure 3. The longitudinal distribution of the sea surface temperature in the equatorial Pacific 
(10"S-10“N) for April 1985 and April 1987. The SST is taken from the NCEP analysis. 

Figure 4. Net downward fluxes at the top of the atmosphere and the surface, and the heating of 
the atmosphere averaged over the period January-August 1998. Units are W m^. 
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